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In this work, we perform parameter tuning with dynamical anisotropic clover lattices using the 
Schrodinger functional and stout-smearing in the fermion field. We find that ^r/^o is relatively 
close to 1 in our parameter search, which allows us to fix <^o in our runs. We proposed to determine 
the gauge and fermion anisotropy in a Schrodinger-background small box using Wilson loop 
ratios and PCAC masses. We demonstrate that these ideas are equivalent to but more efficient 
than the conventional meson dispersion approach. The spatial and temporal clover coefficients 
are fixed to the tree-level tadpole-improved clover values, and we demonstrate that they satisfy 
the nonperturbative condition determined by Schrodinger functional method. 
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1. Introduction 

Lattice quantum chromodynamics (QCD) has successfully calculated many of meson spec- 
trum; however, there remain many challenges for the lattice community to resolve the myriad 
states present in QCD. The case of the nucleon spectrum is one such battlefield. Consider the low- 
est three states in the N spectrum (N, N' (Pu) and N* (Sn)), for example. Many earlier quenched 
lattice QCD calculations ^ find a spectrum inverted with respect to experiment, with N' heavier 
than the opposite-parity state N*. Although the Kentucky group [^] managed to find the correct 
mass ordering around pion mass 300-400 MeV (after taking care of the effects of the quenched 
"ghosts"), no other lattice group has been able to reproduce the experimental ordering using dif- 
ferent approaches. Furthermore, these are just the lowest few states in the N spectrum. There are 
many more states seen in experiment for which lattice calculations could help to identify particle 
properties. 

This situation suggests an urgent need for full-QCD simulations that can resolve some of these 
issues. In order to get better signal for the excited states (especially for the higher-excited nucleon 
spectrum), one needs a lattice with a fine temporal lattice spacing. At the same time, we also 
want to avoid finite-volume effects. Current dynamical lattice gauge ensembles manage to have a 
reasonable lattice box with spatial dimensions about 3 fm, but unfortunately the lattice spacing is 
about 2 GeV, which is not fine enough to allow determination of more than one excited state. One 
solution to this situation would be to generate anisotropic dynamical lattices. 

The anisotropic lattice has been widely adopted in lattice calculations. It was first used to sim- 
ulate heavy-quark physics, such as charm, back in the era when the lattice spacings were too coarse 
to use the relativistic quark action to simulate heavy quarks. Another main application is for calcu- 
lations, such as glueballs[|3|] and multiple excited-state extraction^], where the anisotropic lattice 
has great improvements over isotropic due to the finer lattice time spacing. Previous results on the 
anisotropic lattice include two-flavor anisotropic dynamical simulations done by CP-PACS[|8|] and 
TrinLat collaboration^]. 

In this work, we will use a three-flavor Sheikholeslami-Wohlert (clover) action with stout- 



link smearing (in the spatial direction only) in the Schrodinger functional scheme [|13|]. We can 
determine the gauge anisotropy by looking at the Wilson loop ratios. The ratio of the PCAC-current 
quark mass with background field in space and in time tells us about the fermion anisotropy, which 
has been determined in the past using the meson dispersion relation in large volume. The rest of 
the coefficients are set to tree-level tadpole improved values, with the tadpole factors determined 
from numerical simulation. The clover coefficients are fixed at tree-level tadpole-improved values, 
which later we demonstrate are consistent with nonperturbative ones determined in the Schrodinger 



functional scheme. Our configurations have been generated using the Chroma[]19[] HMC code with 
RHMC for the third flavor and multi-timescale integration. For more details, see Ref. [11]. 

2. Methodology and Setup 
2.1 Action 

We use tree-level tadpole-improved O (a 2 ) -improved Symanzik gauge action: 

4 = TT I A I fo^y + c x m ss \ + % £ [co&st + c A] 1 , (2. l) 

where the £o is the gauge anisotropy. We adopt the clover fermion action 
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a t Q F = —\ u,m Q + v t W t + t L J\W s -- c' sw Y O ts F ts + ^ Y a ssf F ss , 
u t SOT 1 s so s < s , 



(2.2) 



The v t is a redundant parameter and set to 1 to resolve the doubler problem of Wilson-type 
fermions. The v s reflects the ratio between fermion and gauge anisotropy. The Cg' w are the spatial 
and temporal clover coefficients which are set as follows: 

y, , l ( i\ l 

,3 



-sw 



O 5 C 



SW 



1 f 1 

2 v ' + a 



u t uj 



(2.3) 



these selections are discussed in Ref. [|l2|]. The tadpole factors are later set to fixed values taken 
from our early dynamical simulations, which agree amongst themselves within 1-2%. Thus we are 
left to tune the remaining four coefficients: V s , mo, j8. 

2.2 Stout-smeared links 



We will use fhree-dimensionally stout-smeared link variables Q10|] within the fermion action. 
Note that the smearing does not involve the time direction, so the transfer matrix remains physical. 
As with other smearings, one should check the smearing parameter carefully to avoid potentially 
incorrect short-distance physics. In this work, we set p = 0.22 and n p =2 for exploratory study. 

With a nonperturbative determination of the clover coefficients at the target lattice spacing of 
a = 0.1 fm, the scaling violations are about 1% in amy /V a 2 o [JTvI] . These previous scaling studies 



used an isotropic quenched action. On the left-hand side of Figure |T] (where all the points on the 
graph with fixed m n /m p = 0.7), we show the result from CP-PACS's isotropic Wilson scaling be- 
havior; our simulation measurement on the anisotropic lattice shows a more continuum-like scaling. 
When we use stout-link smearing in the Wilson fermion action, we notice dramatic improvement. 
In the clover case (on the right), similar tests are performed, compared with nonperturbative clover 
coefficients. The stout-link smeared clover action point is compatible with the scaling. 

2.3 Schrodinger functional 

The Schrodinger functional[[l3|] has been implemented in lattice QCD since the 1990's. It 
allows us to simulate at lighter pion mass (since the background field lifts zero modes); from 
the PCAC relation we can check how close our csw in the fermion action is compared with the 
nonperturbative value. 

We modify our definition of the quark mass to be 



M(xo,yo) = r(x ) 



r(yo)-r(yo) , . ,, . , ^(jo) - r(y ) , 
s(xq); M (xo,y ) = r {x ) - -jj— -—" 



*'(jo)-*Cvo 

where r and s are obtained from 

r(x ) = 0.25(do + dQ)f A (xo)/f P (x ); s(x ) = 0.5a d dQfp(x )/f P (x ) 



s'(x Q ), (2.4) 



(2.5) 



Ia(p) (with r = 757^ (75) ) is a correlation function of bulk fields (y, y) and boundary fields at t = 
(Jl, 17): 

fo T (t) = (^v^v)E^(y,0rn(v)>/(tf/ - 0- (2-6) 



y.z 



Similar definitions apply to r 1 and s' , but the f^i P \ now involves the other boundary fields at t = 
T,T — 1 0f, 7]') and a sign change. Later in this work, we will calculate the modified quark mass 
to tune our coefficients. 
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Figure 2: a s M s comparison among different spa- 
Figure 1: The scaling behavior of quenched Wilson gauge tial box volumes 
action with Wilson (left) and clover (right) fermion actions. 
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Figure 3: A m i n (2 2 ) measured in simulation with (left) and without (right) a background field. The jc-axis is 
in units of 5 trajectories. 



On the isotropic lattice, one can also determine the csw coefficient in Schrodinger functional 
scheme from the PCAC relation. We require the nonperturbative value of csw to lie at 

AM = M(x ,y )-M'(x Q ,yo)=AM^\ (2.7) 

where AM^ ' is the tree-level mass splitting obtained from a classic background field simulation 
with the same setup of gauge and fermion actions. Previous dynamical works were carried out 



by Alpha with Nf = 2[|14]] using Wilson gauge only, CP-PACS for two-flavor and three-flavor 



calculations with Wilson and Iwasaki gauge actions [J8J, 15J, |lq ]. However, all applications of the 
Schrodinger functional so far have been limited to isotropic lattices. This work is the first to apply 
the Schrodinger functional to dynamical anisotropic lattices. 

In this work, we will implement the Schrodinger functional setup on anisotropic lattices for 
the first time in dynamical simulations. We find that in the three-flavor (anisotropic) clover action 
(with parameter mo = —0.054673, V s = 1) simulation, the lowest eigenvalue of Q^Q is lifted by 
the background field. This helps to reduce the frequency of exceptional eigenvalues, as shown in 
Figure ||. The eigenvalue changes dramatically from 0-6000 trajectories without the background 
field (right figure). A reduction in near-zero eigenvalues gives us better acceptance rate in the 
HMC. 

We also implement the background field not only in the "?" direction (as conventionally used 
in Schrodinger functional) but also the "z" direction. Less is known about putting the background 
field in the "z" direction. (We label the modified masses M s and M t with respect to the direction 
of the boundary field direction and similarly for the mass difference AM S and AM,.) We measure 
M (jco) with various L z , showing results in Figure ^[ When we increase the length in the z direction, 
good signal appears from 12 to 16 but not beyond 24. This is because the background field becomes 
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Figure 5: The measured M s j (in units of a^ 1 ) at each 
V value as a function of mo 



too weak at large L z . When we increase the background field signal, which is proportional to <J> 
at L z = 24, the signal is still not as good as for L z = 16. Similar checks can be done regarding 
the size of L v v . As we increase the value from 8 to 12, the signal shows improvement. The right 
panel of Figure || shows an optimal choice of 12 3 for the spatial volume; this is what we use in the 
remainder of this work. 

2.4 Proposed conditions 

In summary, we implement Schrodinger functional with background fields in two directions: 
"t" and "z", and we measure the quantities M t s and AM SJ . We will determine the gauge anisotropy 
^ R from the ratios of static quark potential (more details in the following section). We will get the 
fermion anisotropy v s by measuring the PCAC mass ratio in finite-volume Schrodinger functional 
scheme; in our coefficient definition, we are looking for a s M s /a t M t = We set the two clover 
coefficients (cg W ) to the stout-smeared tadpole coefficient, where the tadpole factor is numerically 



tuned. We further check the nonperturbative conditions in Eq. 2.7 from our measured the PCAC 
mass difference AM Sjt ; we find that they are consistent. 

3. Numerical Results 

In this work, we fix j8 = 2 for an exploratory study on the tuning of action parameters. 
The natural way to look for anisotropy in the gauge sector is to start with the static quark po- 
tential. We calculate ratios of Wilson loops involving the temporal direction W s t and those without 



R ss (x,y) 



W„(x+l,y) 



'asym. 



; Rst(x,t) 



W st {x+\,t) 



'asym. 



a a s V s (ta,) 



(3.1) 



The finite-volume differences in the ratio R st and R ss are the same if N t = £rN s . Thus, there 
are no finite-volume differences in V s (ya s ) or V s (ta s /E,R) either. Naturally, one should impose 
R ss (x,y) = Rst(x,t) to get the renormalized E,r. The ^r/^o is consistent with 1 within a few percent 
from our parameter searching. (Figure || shows a special case when v v = 1). We can set = 3.5 
and tune only the remaining two parameters. 

We search the remaining two-dimensional space by varying the parameters v s and mo; the 
corresponding PCAC mass (in units of a^ 1 ) is shown in Figure ||. It is interesting to note that when 
V s is around 1 (the classic value), M s as M t among range of different mo values. Figure || shows one 
of our parameters when mo is fixed close to the chiral limit (in our study, it is at —0.07, slightly 
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Figure 6: The measured M SJ in units of a t as a func- 
tion of v at fixed mo = —0.07 



Figure 7: Meson dispersion relation 



below the critical mass). As we can see, at v s = 1 the a t M s is automatically equal to a t M t . This 
suggests that with our setup (stout-smearing, tree-level tadpole coefficients for the fermion action) 
the tunings are consistent with the classic prediction. 

One question is: how does our condition for the fermion anisotropy compare with the conven- 
tional dispersion relation (c 2 = (E 2 — m 2 ) / p 2 ) in a large volume? To demonstrate that the two are 
consistent at 0(a), we pick one of our simulation points, mo = —0.054673 and v v = 1.0 (where 
a s = 0.1 16(3) fm), and we perform a dynamical simulation on a bigger volume, 12 x 128 without 
Schrodinger functional scheme. We measure the pion dispersion relation (shown in Figure [7]), and 
we find c 2 = 1.088(8). This is about the same amount of discrepancy as in the Schrodinger func- 
tional measurement of M s and M t . Therefore, since it uses smaller volumes, probing the condition 
Ms = M t is a more efficient way to tune the fermion anisotropy v s than the dispersion relation and 
with Schrodinger functional, one can work on smaller pion masses (even near chiral limit) without 
too much additional cost. 

The final check is: how good is our initial tadpole-improved Cg w ? In the Schrodinger func- 
tional scheme, such nonperturbative coefficients are determined by requiring that 

AM = M(2r /4, T /4) - M'(2T /4, T /4) = AM Tree ' M=0 (3.2) 

be satisfied. The tree-level AM value is obtained from simulation in a classical background field. 
In the dynamical simulation with parameter V, = 1 and mo = —0.056266, we find M s = M t is 
satisfied. (See Figure. ||] in which the thickness represents the number of the measurements done 
with that specific choice of parameters.) We further check the NP condition: the tree-level values 
of a s AM are -0.00056166 and -0.00028645, and the measured a s AM Sit are -0.000257(424) and 
0.00009(12). Another case at parameters v s = 0.9 and ra = —0.069 gives M s consistent with 

S t 

M t within the errorbar. The NP c sw condition is also satisfied. We found that tadpole-corrected 
tree-level coefficients with the stout-link smearing are consistent with the nonperturbative 0(a)- 
improved coefficients in the three flavor dynamical simulation. 

By imposing the condition of equality on the two PC AC masses (M t , M s ) measured under the 
imposed background fields in two directions: t and "z", we will find the correct parameters v s and 
niQ. A linear interpolation M i / (v s ,mo) = b sf + d Sit v s + e^,mo would work for small variations in the 
parameter space. Some additional runs in the future will help us better determine these coefficients; 
a good interpolating form will help us when tuning the strange quark in the future. 
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4. Conclusion and Outlook 

We demonstrate that the Schrodinger functional combined with stout-link smearing shows 
promise in dynamical runs. Further, we show in Schrodinger functional scheme that stout-link 
smearing and nonperturbatively modified tadpole factors automatically make our 0(a) -improved 
coefficient tuning condition fulfilled. Our proposed finite-box fermionic anisotropy v s tuning 
(using the ratio of the PCAC mass measured with background field in space and time directions) is 
as good as conventional large-box runs but more efficient. 

Our generation of two-flavor anisotropic = 3) Wilson fermion configurations is complete. 
This includes two lattice sizes: L m 1.8,2.6 fm with m n « 400,570 MeV. In the near future, we 
will begin to fine tune the strange quark mass. We will also calculate 0(a) -improved coefficients: 
cva and Zy^ for people who are interested in using these configurations for other physics. 
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